The Evolving ISM in the Milky Way & Nearby Galaxies 



Tracing the [Fell] /[Nell] ratio and its relationship with other ISM 
indicators within star forming dwarf galaxies: a Spitzer IRS 

archival study. 

B. O'Halloran 1 , S. C. Madden 2 , N. P. Abel 3 
boh@physics.gmu.edu 

ABSTRACT 

Archival Spitzer observations of 41 starburst galaxies that span a wide range 
in metallicity reveal for the first time a correlation between the [Fell] /[Nell] 
26.0/12.8 /im ratio and the electron gas density as traced by the 18.7/33.4 /im 
[SIII] ratio, with the [Fell] /[Nell] ratio decreasing with increasing gas density. 
We also find a strong correlation between the gas density and the PAH peak to 
continuum strength. Using shock and photoionization models, we see the driver 
of the observed [Fell] /[Nell] ratios is metallicity. The majority of [Fell] emis- 
sion in low metallicity galaxies may be shock-derived, whilst at high metallicity, 
the [Fell] emission may be instead dominated by contributions from Hn and in 
particular from dense PDR regions. However, the observed [Fell] /[Nell] ratios 
may instead be following a metallicity-abundance relationship, with iron being 
less depleted onto grains in low metallicity galaxies - a result that would have 
profound implications for the use of iron emission lines as unambiguous tracers 
of shocks. 

Subject headings: galaxies: ISM — infrared: galaxies — infrared: ISM — ISM: 
dust, extinction — ISM: structure 
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1. Introduction 



The presence of massive stars within starbursts undoubtedly plays a huge role in de- 
termining the physical conditions within the local interstellar medium (ISM). High-mass 
(M init > 8M ) stars formed within typical starbursts drastically affect the dynamics of the 
surrounding ISM, through not only the release of ionizing photons which destroy molecular 
material, but also via supernovae (SNe) which provide thermal and kinetic energy input 
into the ISM. The effects of photoionization by high-mass stars on the observed dearth of 
PAHs have recently been investigated (Madden et al. (2006); Wu et al. (2006)), while in 
O'Halloran et al. (2006) (hereafter Paper I), we examined a sample of 18 galaxies of varying 
metallicity (from l/50th to super solar) with high star formation rates in order to determine 
whether supernova-driven shocks do indeed play a role in the PAH deficit in low metallicity 
environments. If we consider the ratio of the 26 jum [Fell] line and the 12.8 /im [Nell] line 
as a tracer of the strength of supernova shocks, we found a strong anti-correlation suggest- 
ing that strong supernova-driven shocks are indeed present within low metallicity galaxies. 
Furthermore, the PAH deficit within these objects may indeed be linked to the presence and 
strength of these shocks. However, it has not as yet been conclusively proved that shocks 
are the dominant process behind the PAH deficit. It would therefore be advisable to fur- 
ther explore the nature of the ISM within such star-forming environments using additional 
mid-IR probes in order to expand upon our understanding of the [Fell] /[Nell] ratio, and by 
extension, its relationship with key ISM indicators such as PAH strength and metallicity. 



2. Sample and Data Reduction 

To accomplish this, we have expanded our sample from the 18 objects presented in Paper 
I by including archival IRS observations of 23 additional objects, bringing the sample total 
to 41. These galaxies range in metallicity from extremely low (such as I Zw 18, with Z /Z Q = 
1/50) to super-solar metallicity (> 1 Z®) galaxies such as NGC 7714. None of the galaxies 
in our sample are known to harbour AGNs. We extracted low and high resolution archival 
spectral data from the Short-Low (SL) (5.2 - 14.5 /im), Short-High (SH) (9.9 - 19.6 /im) and 
Long-High (LH) (18.7 - 37.2 /im) modules of the Spitzer Infrared Spectrograph (IRS). The 
datasets were derived from a number of Spitzer Legacy, GO and GTO programs released to 
the Spitzer Data Archive, and consisted of either spectral mapping or staring observations. 
We obtained fluxes for the nuclear positions only from the mapping observations. All the 
staring observations were centered on the galaxy's nucleus. 
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3. Results and Discussion 

In Fig 1, (top left) we plot [FeII]/[NeII] vs the 6.2 /im PAH peak to continuum ratio 
to check if the relationship between the [Fell] /[Nell] ratio and the PAH peak to continuum 
strength first noted in Paper I. Employing a Spearman rank correlation analysis to assess 
the statistical significance of this trend yields a correlation coefficient of r s of -0.705 and 
P s of 2.69 x 10 -6 , confirming the significant anti-correlation between PAH strength and the 
[Fell] /[Nell] ratio as seen in Paper I. We plot the [Fell] /[Nell] ratio versus the metallicity 
in Fig.l (top right), and we again see a strong anti-correlation (r s of -0.777 and P s of 9.36 x 
10~ 7 ) between [Fell] /[Nell] and metallicity as seen in paper I. 

Based upon what we have already seen from Paper I and from the figures, if the 
[Fell] /[Nell] ratio is truly indicative of the strength of supernova- driven shocks within the 
extent of the high resolution slits, one would expect the passage of such shocks to affect con- 
ditions within the local ISM in quite a substantive manner. Shocks, in addition to removing 
dust and PAH from the ISM, should also be adept in removing gas - one would therefore 
expect the propagation of intense SNe-driven shocks into the ISM of star forming regions to 
greatly impact on the density of the gas. In order to determine how the gas density within 
these nuclear star forming regions corresponds with the strength of the supernova-driven 
shocks, we require a reliable mid-IR tracer to probe the gas density. The [SIII] 18.7/33.4 /im 
line ratio provides such a reliable mid-IR tracer of the gas density, as it is ideal for probing 
gas (with high critical densities) in the regions surrounding starbursts, especially as the [SIII] 
lines are starburst dominated. This ratio is sensitive to changes in the density for the 50 < 
n e < 10 4 cm -3 , but is insensitive to changes in temperature. We also plot the logarithm of 
the ratio of the [SIII] fluxes versus the [Fell] /[Nell] flux ratio (bottom left). There is a strik- 
ingly strong trend between the two line ratios, with high [Fell] /[Nell] values corresponding 
to lower [SIII] ratios. Using the Spearman rank correlation analysis, we get (r s ) of -0.891 
between the [SIII] ratio and the [Fell] /[Nell] with a probability of chance correlation (P s ) of 
7.43 x 10~ 8 , indicating a significant anti-correlation. Interestingly, we see a general decrease 
in the metallicity of the object with lower [SIII] ratios, corresponding to higher [Fell] /[Nell]. 

We again plot the logarithm of the ratio of the [SIII] fluxes, but this time against the 
6.2 micron PAH peak/continuum ratio (bottom right Fig. 1). We see a strong trend, but 
this time with the PAH strength increasing with log [SIII] - the PAH strength is increasing 
with higher gas density. Again running a Spearman rank correlation analysis, we get y s of 
0.750 and P s of 9.64 x 10~ 7 , confirming a significant correlation between log [SIII] and PAH 
strength. Breaking down the points in metallicity bins, we again see that a general trend 
exists with PAH strength increasing with metallicity and increasing values of log [SIII]. 

Ruling out extinction and aperture effects as drivers for the observed correlations, we can 
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be confident that the correlations seen in Fig. 1 are indeed true physical relationships between 
the strength of the [Fell] /[Nell] ratio and a number of ISM indicators, with the primary 
culprit for the observed relationships being the passage of SNe-driven shocks. However, 
while we have focussed up to now with a supernova-derived shock origin for the behaviour of 
the [Fell] /[Nell] ratio, it may not be the only explanation. For example, Izotov et al. (2006) 
note sign of strong depletion of iron onto dust grains, and gradual destruction of those grains 
on a time scale of a few Myr, based on a survey of metal-poor galaxies from the 3rd release 
of the SDSS. Such a process could undoubtedly drastically affect the nature of Fe emission, 
and by extension the behaviour of the [Fell] /[Nell] ratio, within our sample - we may instead 
be probing an abundance-driven relationship and the other observed relationships presented 
here would be purely coincidental to the [FeII]/[NeII] ratio. In order to explore if this indeed 
is the true cause or if shocks alone can explain the observed relationships, we used shock and 
standard Hn - PDR models in an effort to model the observed [FeII]/[NeII] relationship for 
a wide variety of environments. We used both the MAPPINGS III photoionization/shock 
code and the Cloudy photoionization code in order to determine the relative proportions 
to the [Fell] emission from shocks and Hn/PDRs, and by extension, the driving process 
behind the behaviour of the observed [Fell] /[Nell] line ratio. We found that the relative 
Hn contribution dominates over the PDR contribution at low densities (and metallicities), 
with the PDR contribution becoming more dominant as one moves to higher densities (~ 
10 3 cm -3 ) and metallicities. For very dense Hn regions where the PDRs are irradiated 
by intense FUV, and thus PDR-derived [Fell] emission is dominant. Correspondingly at 
lower densities (~10 2 cm -3 ), the Hn regions are larger, the PDRs lie further from the 
stars, and the resultant lowered FUV fluxes and densities do not excite [Fell] in the PDR, 
and thus the Hn region dominates the Hn/PDR contribution to [Fell] production for such 
environments. From the MAPPINGS and Cloudy simulations and the comparisons with the 
IRS and optical data, it would seem that metallicity is the primary driver for the behaviour 
of the [FeII]/[NeII] ratio, as metallicity determines the origin (shocks vs Hn and/or PDRs) 
and strength of the [Fell] emission. We can reproduce the observed [Fell] /[Nell] dependence 
on Z if we assume the Izotov variation of Fe/Ne with Z. However, we do note that the Hn 
and PDR contribution to the [Fell] emission is surpassed by the contribution from shocks 
at the lowest metallicities and densities, based on comparison with the MAPPINGS output 
and the observed data - perhaps only ~10% of the [Fell] emission comes from H II and PDRs 
with that contribution rising to ~80-90% as one moves to higher metallicity. As we move to 
roughly solar metallicity, it would seem likely that the dominant mechanism for production 
of [Fell] (and thus the behaviour of the [Fell] /[Nell] ratio at high metallicities) is from Hn 
and PDRs in combination rather than from SNe shocks. 
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Fig. 1. — (Top left)Plot of the [FeII]/[NeII] ratio versus the PAH 6.2 micron peak to con- 
tinuum ratio for the extended sample. (Top right) Plot of the [FeII]/[NeII] ratio versus the 
metallicity for the extended sample. (Bottom left) Plot of the 18.7/33.4 /mi [SIII] ratio as 
a function of the [Fell] /[Nell] ratio for the extended sample. (Bottom right) Plot of the 
18.7/33.4 /ini [SIII] ratio as a function of the 6.2 /mi PAH peak to continuum strength for 
the extended sample. 
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4. Summary 

Using archival Spitzer observations of 41 starburst galaxies that span a wide range in 
metallicity, we found a correlation between the ratio of emission line fluxes of [Fell] at 26 
/im and [Nell] at 12.8 /xm and the electron gas density as traced by the 18.7/33.4 fim [Sill] 
ratio, with the [Fell] /[Nell] flux ratio decreasing with increasing gas density. We also find 
a strong correlation between the gas density and the PAH peak to continuum strength. 
The correlation of the [Fell] /[Nell] ratio and the PAH peak to continuum strength found 
in paper I was confirmed for a larger sample. Using shock and photoionization models, we 
see that metallicity is the primary driver for the observed behaviour of the [Fell] /[Nell] 
ratio. It may very well be that the majority of [Fell] emission at low metallicity may be 
shock-derived, whilst at high metallicity, the [Fell] emission is dominated by contributions 
from Hll and in particular, from PDR regions, and that at higher metallicity shocks may 
not play as significant a role in removing gas, PAH and dust from the ISM, unlike in low 
metallicity systems. However, the observed [Fell] /[Nell] emission may instead be following 
a metallicity- abundance relationship, with the iron being less depleted in low metallicity 
galaxies, a result that would have profound implications for the use of Fe emission lines as 
unambiguous tracers of shocks. 
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